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On  Singular  Semilinear  Elliptic  Equations 

Aihua  W.  Shaker 
Department  of  Mathematics 
Naval  Postgraduate  School 
Monterey,  CA03O43 

May  8,  1991 


Abstract:  -  For  the  semilinear  elliptic  equation  Au  +  p(x)u~"'  =  0,  x  €  71", 
n  >  3,  7  >  0,  we  show  via  the  barrier  method  the  existence  of  a  positive 
entire  solution  behaving  like  |a:|2""  near  oo. 

1  Introduction 

We  study  the  singular  semilinear  elliptic  equation 
(1)  Au  +  p{x)u""1  =  0 

in  72".  This  type  of  equation  arises  in  the  boundary  layer  theory  of  viscous 
fluids  [3,4].  From  the  results  of  Fulks  and  Maybee  [8],  Crandall,  Rabinowith, 
and  Tartar  [5],  Gomes  [9],  and  recently  Lazer  and  McKenna  [14],  it  follows 
that  (1)  has  a  unique  classical  solution  within  a  bounded  domain  fl,  where 
p(x)  is  a  sufficiently  regular  function  which  is  positive  on  Cl.  Kusano  and 
Swanson  [12]  gave  the  existence  proof  on  exterior  domains.  As  for  the  ex¬ 
istence  of  entire  solutions,  not  much  is  known.  Edelson  [7],  Kusano  and 
Swanson  [13]  have  been  able  to  show  the  existence  of  entire  solutions  of  (1) 
with  7  €  (0,1),  and  p(x)  sufficiently  regular.  In  this  paper  we  show  via  the 
upper  and  lower  solution  method,  which  is  also  referred  to  as  the  barrier 
method,  that  (1)  has  a  bounded  positive  entire  solution  vanishing  at  oo  in 
Rn  for  n  >  3  and  all  7  >  0. 

The  author  learned  after  this  paper  was  finished  that  a  similar  result  was* 
given  earlier  by  R.  Dalmasso  [6],  but  by  a  different  approach. 
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2  Preliminaries 


We  first  state  the  theorem  by  Kusano  and  Swanson  [13]  for  the  case  0  <  7  < 
1. 

Lemma  1.  Equation  (1)  has  an  entire  bounded  positive  solution  t i(z)  in  Rn 
for  n  >  3  ,and  |x|n_2u(x)  is  bounded  and  bounded  away  from  zero  near  00  if 
p(x)  satisfies  the  following  conditions: 

(HI)  p(x)  e  Cfoc{Rn),  n  >  3,  p(x)  >  0,  x  e  Rn  \  {0}, 

(H2)  3C  >  0,  such  that  C<f>( |x|)  <  p(x)  <  <£(|x]),  ip(x )  =  max|x|=(p(x),  0  < 
t  <  00, 

(H3)  /r  tn-i+^n~2^(t)dt  <  00. 

The  term  ’’entire”  has  often  been  used  for  solutions  of  equation  (1)  in 
Rn.  To  avoid  confusion  with  the  traditional  definition  for  entire  functions, 
we  use  the  term  ”C2+a-entire”.  A  C2+a-entire  solution  of  (1)  is  defined  to 
be  a  function  u(x)  €  C?+a(Rn)  that  satisfies  (1)  pointwise  in  Rn. 

The  method  that  we  shall  be  using  heavily  in  our  proof  is  the  so-called 
barrier  method,  or  upper-lower  solution  method. 

We  consider  the  elliptic  boundary  value  problem 

J  Lu  +  f(x,u)  =  0  in  D 

|  Bu  —  adu/du  +  bu  —  g  on  OD 

where  D  is  a  smoothly  bounded  domain  in  i?A  and  1/  =  (// 1,  •••.*'„)  is  a 
smoothly  varying  outward  normal  vector  field  on  dD  which  is  of  class  C2+a, 
while  a  and  b  are  positive  constants.  We  also  assume  that  /  £  Ca  and  that 
g  has  an  extension  g  to  the  interior  of  D  such  that  g  €  C2+a. 

An  upper  solution  to  the  above  problem  is  a  function  <p  satisfying 

f  Lip  +  /(x,  <p)  <  0  in  D 
(  B<P  >  g  on  dD. 

A  lower  solution  to  the  above  problem  is  a  function  ip>  satisfying 

f  Lip  -f  f(x,ip)  >0  in  D 
\  Bxp  <  g  on  3D. 

Wc  assume  that  dD,f,g,  and  the  coefficients  of  L  are  smooth  in  what 
follows 


Lemma  2.  (Theorem  2.3.1  of  [16])  Let  ep  be  an  upper  solution  and  ip  a 
lower  solution  with  ip  <  <p  on  D .  Then  there  exists  a  solution  u  to  the  above 
boundary  value  problem  with  ip  <  u  <  <p. 

We  consider  the  following  example: 

f  u" ’  +  Xu  —  u3  =  0  x  6  (0, 7r) 

\  u  =  0  x  =  0,  7T 

By  the  above  theorem,  if  A  >  1,  then  the  problem  has  at  least  three  solutions. 

Actually,  u  =  esinrr  with  e  small  is  a  lower  solution,  and  u  =  Rx 1/2  with 
R  large  is  an  upper  solution.  Therefore  there  exists  a  solution  u  such  that 
u  <  u  <  u  in  (0, 7r ).  Clearly  —  u  and  0  are  also  solutions  to  this  problem. 

The  following  lemma  on  the  barrier  method  for  D  —  Rn  is  due  to  Ni 
[15]  in  1982.  A  special  case  was  proved  earlier  by  Ako  and  Kusano  [l]  in 
1964.  The  proof  is  standard.  Using  the  well  known  result  on  the  upper-iower 
solution  approach  in  bounded  regions  (see  Sattinger  [16]),  we  first  solve  the 
equation 

Lu  +  F{x,  u )  =  0 

on  Br.  Then  by  letting,  R  —>  oc,  we  obtain  a  solution  on  Rn  by  a  diagonal 
process. 


Lemma  3.  Let  >  u2  in  Rn  be  such  that 


Lu\  -f  / (x,  Uj)  <  0 
Lu2  +  f{x,  u2)  >  0 


where  f  is  locally  Holder  continuous  in  ( x,u )  and  locally  Lipschitz  in  u ,  and 
L  is  an  elliptic  operator  of  second  order.  Then  there  exists  a  solution  u  of 
Lu  +  f(x,  u)  =  0  with  u ,  >  u  >  u2. 


□  □ 


r 'he  difficulty  in  constructing  the  proof  is  to  find  an  appropriate  upper 
solution  c  equation  (1).  In  order  to  use  the  barrier  method  we  first  study 
the  nonsi.  ^ular  equation 

A  u  +  p(:r)[£  +  =  0. 


For  each  fixed  7  there  corresponds  a  solution  Letting  7  — >  00,  we 

show  that  the  limiting  function  is  the  desired  solution. 

Proof:  By  Lemma  1,  for  7  =  71  £  (0,1),  equation  (1)  has  a  C2+a-entire 
positive  solution  ^(2)  in  Rn,  n  >  3,  vanishing  at  00  at  the  rate  r2-n.  We 
claim  that  u  =  cu\  is  an  upper  solution  of  the  equation  (1)  for  7  >  1,  where 


9  < 


1  +  7i 

1  +  7 


<  1, 


(4)  A  11(2)  +  p(i)[u(2)  +  6]  1  =  0,  2  £  Rn. 


u  =  0  is  a  lower  solution  of  (4).  Since  u  =  cu ’  >  0,  u  >  u  in  Rn .  By  Lemma 
2,  (4)  has  a  solution  u  such  that  u  <  u  <  u. 


For  8  <  8,  u  is  a  lower  solution  of  (4)  with  8  =  8.  Lemma  2  then  implies 
that  (4)  has  a  solution  it  for  8  =  8  such  that  u  <  u  <  u. 

Let  {<5n}^  be  a  sequence  of  strictly  decreasing  positive  numbers,  and  let 
u„(t)  be  a  smooth  positive  solution  of  (4)  when  8  =  8n.  From  the  construc¬ 
tion  of  our  lower  solutions,  it  is  clear  that  un(x)  >  u„_i(x)  for  all  n.  So 
lim^oo  un(x)  =  u(x)  exists  for  all  x  £  Rn  and 

(5)  u  <  u  <  u 


for  x  £  Rn. 

We  can  now  assert  that  u  £  C2+a(Rn)  and  that 

(6)  An  +  p(a,)u-"’ =  0 

for  x  £  Rn.  This  follows  from  more  or  less  standard  arguments. 

Let  x0  £  Rn  and  r  >  0.  We  consider  the  ball  of  radius  r  centered  at  x0 , 
B(x0,r)  in  Rn.  Let  ^  be  a  C°°  function  which  is  equal  to  1  on  B(x0,r/ 2) 
and  equal  to  0  off  B(x0,r).  We  have 

A(^un)  =  2VT  •  Vun  -f  pn 

for  n  >  1.  where  pn  is  a  term  whose  L 1X1  norm  is  bounded  independently  of 
n.  Therefore  for  n  >  1  we  have 

A,  d{^un) 

tyunA(tyun)  =  V  bn]~ - -  +  qn, 

%  dxi 

where  bnj,j  =  1  and  qn  are  terms  bounded  independently  of  n  for 

it  >  1.  Integrating  the  above  equation,  we  have  that  there  exist  constants 
ci  >  0  and  c2  >  0  independent  of  n  such  that 

/  |Vun|2cLr  <  Ci(  /  |Vu„|2d.r)1/2  +  c2. 

JB(x0,r)  JB{i0,T ) 

From  this,  it  follows  that  the  L2(B(x0. 7-))-norm  of  | V 'I' u „ !  is  bounded  inde¬ 
pendently  of  n.  Hence,  the  L2(B(x0,r / 2))-norm  of  |Vn„|  is  bounded  inde- 
pendently  of  n.  Let  be  a  6'°°  function  which  is  equal  to  1  on  B(x0.r/ 4) 
and  equal  to  0  iff  B(x0,r/ 2).  We  have  for  n  >  1, 


A('LiUn)  =  2V'I'i  •  V un  +pln. 
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where  pln  is  a  term  w'hose  L°°(B(x0.r /'. 2))-  norm  is  bounded  independently 
of  n.  From  standard  elliptic  theory,  the  W2'2(B(x0,  r/2))-norm  of  $1^  is 
also  bounded  independently  of  n  and  hence,  the  W2'2(B(x0,  r/4))-norm  of 
un  is  bounded  independently  of  n.  Since  the  W1,2[B(x0,  r/4))-norm  of  the 
components  of  Vun  are  bounded  independently  of  n,  it  follows  from  the 
Sobolev  embedding  theorem  that  if  q  =  2n/(n  —  2)  >  2  for  n  >  2  and  in 
addition  if  q  >  2  is  arbitrary  for  n  <  2,  then  the  Lq(B(x0,  r/4))-norm  of  |itn| 
is  bounded  independently  of  n.  Let  $2  be  a  C°°  function  which  is  equal  to  1 
on  B(x0,r/8)  and  equal  to  0  iff  B(x0,r/4).  We  have  for  n  >  1, 

&{^2un)  —  2V'I,2  •  Vun  +  p2n, 

where  p2n  is  a  term  whose  L°°(B(x0,r/ 4))-  norm  is  bounded  independently  of 
n.  Since  the  right  hand  side  of  the  above  equation  is  bounded  in  Lq(B(x0,  r/4)) 
independently  of  n,  the  W2,<J(B(a-0,  r/4))-norm  of  is  also  bounded  inde¬ 
pendently  ol  n.  Hence,  the  W2,q(B(x0,  r/8))-norm  of  un  is  bounded  indepen¬ 
dently  of  n.  Continuing  this  line  of  reasoning,  after  a  finite  number  of  steps, 
we  find  a  number  rq  >  0  and  q\  >  n/(l  —  a)  such  that  the  W2,91  (B(x0,  rx  ))- 
norm  of  u„  is  bounded  independently  of  n.  Hence,  there  is  a  subsequence 
of  {un}f,  which  we  may  assume  is  the  sequence  itself,  which  converges  in 
C1+Q(B(x0,  )).  If  6  is  a  C°°  function  which  is  equal  to  1  on  (B{x0,r1/ 2)) 

and  0  off  B(x0,r  1),  then 

A(6un)  =  2 V6  ■  Vun  +  pn ,  where  pu  =  6Aun  +  unA6. 

The  right-hand  side  of  the  above  equation  converges  in  Ca{B[x0^r-i)).  Hence 
by  Schauder  theory,  {0u„}~  converges  in  C2+a(B(x0,  ))  and  thus  {un}f 

converges  in  C2+o(B(x0,  /2)).  Since  x0  was  arbitrary,  this  shows  that  u  G 

C2+Q{Rn).  Clearly  (6)  holds. 

4  Some  remarks 
Remark  1: 

For  n  =  1,  the  properties  of  positive  solutions  of  equation  (1)  have  been 
studied  by  Taliaferro  [17],  and  Gatica  [10].  For  ri  =  2.  no  entire  positive 
solution  of  equation  (1)  exists  regardless  of  its  asymptotic  behavior  at  00 
(see  [13]). 

Remark  2: 
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It  is  observed  by  Callegari,  Friedman  and  Nachman  [2],  [3,4]  that  if  the 
partial  differential  equations  describing  the  boundary  layer  behind  a  rarefac¬ 
tion  or  shock  wave  (with  viscosity  proportional  to  the  temperature)  traveling 
down,  and  perpendicular  to,  a  flat  plate  are  written  in  terms  of  a  stream  func¬ 
tion  and  a  similarity  variable  the  following  Blasius-type  equation  emerges 
[18]. 

f"'(v)  +  Mf'Xv)  =  o. 

where 

/( 0)  =  0,  /'(0)  =  A\  /'( oc)  =  l. 

Here,  0  <  K  <  1,  for  rarefaction  waves  and,  1  <  K  <  6,  for  shock  waves. 
(I\  =0  corresponds  to  the  classical  Blasius  problem.)  Adopting  the  Crocco 
variables 

2-  =  fXvh  9  =  f'Xvi) 

results  in  the  system 

99  +  a-  =  0, 

g\ A')  =  0,  flr(l)  =  o, 

which  falls  into  the  class  of  equation  discussed  in  this  paper. 
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